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Superconducting double transition in PrOs4Sbi2 is investigated by analyzing the anisotropy 
of the upper critical field Hc2 in the ab-plane, and the possible pairing state is discussed. When 
mixing due to gradient coupling is active, the twofold-symmetric component is necessarily 
induced in the fourfold symmetric phase, leading to the twofold oscillation of Hc2, contrary to 
the experimental result. To avoid the mixing effect, the weak spin-orbit coupling triplet pairing 
state is considered as a likely pairing function, where time-reversal symmetry is broken. 
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The newly discovered superconductor PrOs4Sbi2 is a 
heavy-fermion compound with filled skutterudite struc- 
ture.^^ It is considered that quadrupole fluctuation is 
related to the pairing of superconductivity, suggesting 
a new pairing mechanism. In order to study the pair- 
ing mechanism, it is important to identify the pairing 
symmetry. The experiment on the nuclear spin-lattice 
relaxation rate suggested that PrOs4Sbi2 has full gap or 
point nodes, excluding line nodes, at a zero field. The 
unique phenomenon in the study of the pairing symme- 
try is the superconducting double transition. The spe- 
cific heat jump at the transition temperature is broad, 
suggesting a double transition. By the thermal trans- 
port measurement in the magnetic field (H) rotated in 
the a6-plane of the crystal axes, the phase diagram of 
the double transition was revealed. ^-^ This experiment 
also suggested that a fourfold-symmetric pairing function 
around the c-axis in the high-field H-phase is changed to 
a twofold-symmetric one in the low-field L-phase at the 
second transition, when magnetic field or temperature is 
decreased. Thus, the H-T phase diagram is divided into 
two regions, H- and L-phases. 

To explain the double transition, we need multiple 
components for the pairing functions, i.e., an unconven- 
tional pairing other than the conventional s-wave pair- 
ing. The fourfold-twofold transition in PrOs4Sbi2 can 
be explained as follows: After the fourfold-symmetric 
pairing component appears at the first transition, the 
second component appears at the second transition, and 
the combination of these components gives the twofold- 
symmetric gap structure. For PrOs4Sbi2, the scenario 
of "anisotropic-s"-|-id-wave pairing^^ with point nodes^^ 
was proposed. This can explain the double transition at 
zero field under the cubic symmetry. However, in this 
scenario, it is difficult to reproduce the fourfold-twofold 
transition in magnetic fields, as will be discussed later. 

In the vortex state at a finite magnetic field, we have 
to consider the mixing of the multiple components of 
the pairing function. In the multicomponent Ginzburg- 
Landau (GL) theory, this mixing comes from the gradi- 
ent coupling between the different pairing components. 
This mixing effect has been well studied in the d-wave 
pairing case for high-Tc superconductors, ''^^'^^ and in the 
triplet pairing case for Sr2Ru04 and UPts.^''"^^^ In the 



magnetic field, if the mixing effect is active, a twofold 
component is induced even above the fourfold-twofold 
transition, indicating that the resulting superconduct- 
ing gap shows a twofold-symmetric character and that 
the phase transition between the H- and L-phases is 
smeared. The easiest way to clarify the effect of the 
twofold-symmetric character is to study Hc2 anisotropy 
in the magnetic field rotated in the afo-plane. If Hc2 
has twofold anisotropy, all quantities should show the 
twofold anisotropy near Hc2- Since twofold anisotropy 
is not observed in the experiment, we have to consider 
the condition that it does not occur. This gives impor- 
tant information for identifying the pairing function for 
PrOs4Sbi2. Also in UPta, a theoretical study to explain 
the phase diagram of the double transition in the mag- 
netic field strongly suggested triplet pairing. 

The purpose of this letter is to study Hc2 anisotropy 
for the fourfold-twofold transition case, and to estimate 
the effect of the mixing between fourfold and twofold 
pairing components. By obtaining an analytical result 
for Hc2 near Tc in the microscopic theory, we show gen- 
erally that Hc2 shows twofold oscillation, contrary to the 
experimental result if gradient coupling terms are in- 
cluded. Therefore, gradient coupling terms should be 
absent in order to assure the fourfold behavior above the 
second transition. Our calculation is performed using 
the general pairing function forms for both singlet and 
triplet pairings, so that our analysis can be applied to 
any pairing function cases. In the last part, we discuss 
the possible pairing function for PrOs4Sbi2. 

For the multicomponent superconductors, the pair po- 
tential is decomposed to 



A(r,fc)Q/5 



(1) 



with the order parameter rjm{r) for the m-th pairing 
component. Here, r is the center of mass coordinate 
of the Cooper pair. The relative momentum k of the 
pair is mapped on the Fermi surface. The gap function 
is a 2 X 2 matrix of a and (3 (=t, i), given by 







0m(fe) 



(2) 
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for singlet pairing, and 



-dr, 



dm.: 



(fe) 



(3) 



for triplet pairing. The pairing function satisfies the nor- 
malization and orthogonahzation condition {<j)^(j)m')k = 
dm,m', where fe^ m' ^ t r((^]:^(^^/)/2 = (t)^^m' for sin- 
glet pairing and 



^m.x^rn' ,x + d. 



+ 



djjj • dm' for triplet pairing. {■ ■ ■)k rneans 
the average of fe on the Fermi surface. Similarly, the 
pairing interaction is decomposed to 



V{k',k)a'0',a0 = Vm4>Tn{k')a'p'$ln{k)ai3. 



(4) 



The transition temperature Tc^m for the m-th component 
is given by the relation (V^A^f)"^ = ln(2a;c7/7rrc,m) 
with the Eulcr constant 7 and the cutoff frequency cjc. 
A''f is the averaged density of states on the Fermi surface. 

We derive the gap equation along the Hc2 line by the 
quasiclassical method. ^'^^ The 2x2 matrix quasiclassi- 
cal Green's functions / and g are determined using the 
quasiclassical Eilenberger equation 

i2L0n + iv ■ q)f = (gA + Ag), 

{2u;n + ivq*).f = {gA^ + A^g), (5) 

where g^ + ff = 1 and g^ + ff = 1 with a unit 
matrix 1. In eq. (5), v is the Fermi velocity, and 
q = ~\{d/dr) + {2tt / (pi^) A{r) with a vector potential 
A{r) = {—Hy, 0, 0) in Landau gauge and the flux quan- 
tum 00- The selfconsistent condition is given by 

A(r, fe')a'/3' = — 

a/3 \aj„\<u!c 

Xf{r,k,iWn)al3)k- (6) 

After eq. (5) is linearized about A along Hc2, 
we expand / in powers of v ■ q,^^^ giving faff = 
E~o{-isgnK)/2p|a;„|-(^+i)(t; • q)^ A„^. By substi- 
tuting this in eq. (6), we obtain the gap equation 



Vm' = NpVm' A2n ■ qf\ 



■i)h'nm, (7) 



n=0 



with A2n = 2(-l)"(27rr)-2"(i _ 2-(2"+i))((2n + 1) for 
n > 0, and ln(2cjc7/7rT) = {VxN^)-^ +\n{T^^i/T). 

We assume that the first component 771 is the dominant 
one with the highest transition temperature Tc = T^^i, 
and that other components (m > 2) appear at the 
lower temperature Tc_m < i.e., double transition 
occurs at zero field. We solve eq. (7) in perturba- 
tion of ln(rc/T), and obtain Fc2 near T^J'^^ By in- 
troducing annihilation and creation operators for the 
Landau level as a = —e^^^^{qx — ^^^^Qy) ^md — 
—€~^^^{qx + iiy^^Qy), respectively, we can write v ■ q = 
—v-F^{vj[^a + u_a^), where e — AirH/cfjoi', v± — {vx ± 
ii'Vy)/2, V = v/v-p with the averaged Fermi velocity 
vp- Since we assume that z-axis is parallel to the ap- 
plied field, Qz = 0. The anisotropy ratio u is defined 



as 



= {{Vl\mk/{v'y\cl^l\')ky/' with |<^i|2 ^ 



1, so 



that {Vj.\<pi\'^)k = 0. In dimensionless form, eq. (7) is 
written as 

= ryi - 2C^e ^((^-^a + V-a)f'^^^)r]m 

m 

+2^{eef ^((t)+a + v-a)f'^)nm + • • • , (8) 

m 

= -?7m' - 27C„/^2£^^(^^^^^_^t)2^*^^^^^^^ 

m 

+27'c,„'(e'e)' Y.{{v+a + V-a^f'^;^)r]^ + ■■■ (9) 

m 

for m' > 2, where 7 = a\n{TjT), a = 2|A4|/|^2|^ = 
0.908, e = (|^2h^|/21n(Te/T))i/2 and l/c^ = 
ah).(T /Tc^m)- Equations (8) and (9) include the gra- 
dient coupling terms between different components and 
the nonlocal correction terms in the GL theory. In terms 
of the Landau level function ipn{r) = (r7,!)~^/^(a^)"^o('') 
for the n-th Landau level, the order parameters are de- 
composed to (r) = J2'^=oVm,n'(pn{r). By expanding 
r)m,n and e in powers of 7, we can solve eqs. (8) and (9) 
in perturbation. As a result, we obtain 



Hc2 = 



4n^ 



^1/2 



X {1 + 7(ei + J2 + 0(7^)} 



^>2 



with 



ei = 



2(fi+^_|0iP)|' 



(10) 



(11) 



{V+V. 



(m > 2) (12) 



Along Hc2, the order parameters are given by 

771 (r) = 77i,o[V'o(t') + 'y{fii,2tp2{r) -|- 77i,4V'4(r)} + 0{j^)], 

(13) 

Tlmir) = CmVlAliVrnfii^oir) + fim,2lp2{r)} + 0(7^)] (14) 

for 771 > 2. The factors rim,n are obtained as in eqs. 
(2.12)-(2.16) of ref. 7. 

We consider the case when a magnetic field is applied 
parallel to the a6-plane of the crystal coordinate. In 
this case, the transformation of the Fermi velocity from 
{va,Vb,Vc) in the crystal coordinate to {vx,Vy,Vz) in the 
vortex coordinate is given as {vx,Vy,Vz) = {—VaSinO + 
Vb cos 9, Vc,Va cos 9 + if, sin 9), where 9 is the rotating 
angle of the magnetic field in the a6-plane. We as- 
sume that the superconducting gap for the dom- 
inant component has a fourfold-symmetry around the c- 
axis, so that (Wal«^iP)fc = {vb\<Pi\'^)k and {vaVb\(pi\'^) k = 



0. Th erefore, we obtain {vl\ 4>i\'^) k = JvlH'i\'^)k, 
{vl\cl^i\^= {vl\mk, = {iM')k/{v'c\4>i\')k and 
{v+v-\(l)i\'^)l = {vl\(j)i\'^)l/4. As a result, 

3 



ei = 



^{€\H%'{Mivi+'^'vir\^i?)k 
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3«^«e')l<^iP)fe(l-cos40)}, (15) 



+2 



-{{K - %)4>l4>m)k COS 29 - 2{vaVb4>l(tim)k sin 26^ 

±\Av{^{VbVc(j)l(j)m)kCOSe - {VaVc(l>l(l>rn)kSm.e^\^ },(16) 



(^M'')k, {vlvl\mk = 

<l>l?)k = (WaUf|0l|2)fe = 



where (;Oa|(^iP)fe 

({igvg|0 ip)fc and 

{VaVbVc\4>l\^)k = 0. 

Prom eqs. (10), (15) and (16), we consider the 6 depen- 
dence of Hc2 when the magnetic field is rotated within 
the a6-plane. The term ei comes from the nonlocal cor- 
rection term in the GL theory, and gives the fourfold 
anisotropy of Hc2 under field rotation because of cos4^? 
in eq. (15), reflecting fourfold-symmetry of the dominant 
pairing function (^i(fc). This fourfold Hc2 anisotropy 
was studied in the d-wavc; pairing case in high-Tc su- 
perconductors and the anisotropic s-wave pairing case 
in borocarbide superconductors/' ^^'^^^ The term Cm 
comes from the mixing of 0i and (f>m due to gradient 
coupling. By this mixing, T]m{r) (m > 2) is induced 
even at T > Tc^m in the vortex state at a finite field, 
as shown in eq. (14), while it does not appear at a zero 
field. When the pairing function (pm [m > 2) has twofold 
anisotropy to introduce the fourfold-twofold transition 
at a zero field, gives twofold anisotropy in the 9 de- 
pendence of Hc2- As an example, we consider the case 
of the "anisotropic s"-|-id-wave singlet pairing scenario 
proposed for PrOs4Sbi2,^^ whose pairing functions are 
given as 




(2^6 — H~ ^a)> 



(17) 
(18) 
(19) 



with k = k/k-p. ' Numerical factors come from the 
normalization condition {\4>m\'^)k = 1- We assume an 
isotropic sphere Fermi surface. If r]3 or 772 appears at the 
second transition into the L-phasc, the fourfold-twofold 
transition is explained at a zero field within point-node 
pairing functions. However, we obtain 



1197 
T15" 
169 



^(l-cos4^). 



(20) 



€2 = 7;7^7;;:7;{(1-3C0S26I)^-|-18(1-|-C0S2(9)^}, (21) 

-{(H-cos26')2 -F2(3-cos26')2}, (22) 



38720 ' 
507 



es = 



38720 ' 



since ((t)^ — vl)4'icf)m)k 7^ 0. Because of cos 261 in £2 
and 63, ffc2 shows twofold anisotropy, indicating that 



the H-phasc becomes twofold-symmetric in the vortex 
state by the induced 772 and 773 components. This is con- 
trary to the experimental result.^^ Approaching T^m, 
the twofold-symmetric contribution of Cm increases be- 
cause of the factor Cm in eq. (10). 

To suppress the twofold oscillation in Hc2, it is neces- 
sary to satisfy the relations 



((*a - '"b)'t>l(t>m)k = {VaVb(f>l(f>m)k 



= {VaVc(l>*i(pm)k = {VbVc(pt<l>m)k = (23) 

in eq. (16), indicating that all gradient coupling terms in 

the GL theory vanish. If we consider the correction terms 
with a higher order of 7, other terms with higher order 
combination of v appear as constraints to be zero. The 
relations in eq. (23) give strong restriction for identifying 
the pairing function. For singlet pairing, it is difficult 
to satisfy all relations in eq. (23), when is fourfold- 
symmetric and is twofold-symmetric. For the triplet 
pairing case, if di ■ d*^ ^ 0, the difficulty in satisfying eq. 
(23) is the same as in the singlet pairing case. However, 
all relations in eq. (23) are fulfilled in the case of di • 
'Kn — rf- vectors of the triplet pairing functions are 

orthogonal to each other. 

We now disciiss the pairing function appropriate for 
PrOs4Sbi2. The possible components of the pairing func- 
tions are classified under the crystal group consideration. 
For the present CTibic symmetry (O/i), complete classi- 
fication was carried out by Volovik and Gor'kov^^^ in 
the strong spin-orbit coupling case, while Ozaki et al. 
enumerated possible p-wave pairing^'"'^ and higher angu- 
lar momentiim states (/-wave, • • • )^®) in the weak spin- 
orbit coupling c;ase. These constitute a list of the possible 
states allowed under Oh for both singlet and triplet cases. 
The present task is to select the possible states from the 
list of classifications under cubic symmetry in the light 
of various theoretical and experimental constraints. The 
constraints we require on the H-phase are (i) no gradient 
coupling by satisfying eq. (23), (ii) point nodes on the a, 
b and c axes,'''^^ and (iii) 0?i-symmetric pairing function. 

Let us start with the singlet category. There are sev- 
eral states with point nodes. For example, (/>i(fc) = 
kl +ekl +e'^k1 [e = e='='2'r/3) ji^s point nodes, but along 
(1,1,1) and equivalent directions. It breaks constraint 
(ii). If we consider ad hoc states accidentally combined 
with higher angular momentum states, we can find other 
point node functions satisfying (ii), such as the s-|-(?-wave 
state given in eq. (17). However, for these singlet states, 
the relations in eq. (23) are not satisfied when we con- 
sider the twofold-symmetric second component. There- 
fore, generally, there is no appropriate state among the 
singlet category free from constraint (i). 

As for the strong spin-orbit triplet category classified 
by Volovik and Gor'kov,^*^'^^) the pairing function with 
point nodes is di{k) — xdi^x+ydi^y+zdi^z = xka+y£kb+ 
ze^kc in the p-wave pairing. This is a non- unitary state, 
where one branch of the energy gap jdip ± jdi x is a 
full gap, and the other has 8 point nodes along (1, 1, 1) 
and equivalent directions, contrary to constraint (ii). In 
/-wave pairing, we have di(fc) = xkaik"^ — kf) + ykbik"^ — 
kl) + zkc{kl-kl) and di(fe) = xka{kl-kl)+yekb{kl- 
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fc^) + ze^kcik^ — fc^). These have 14 point nodes, i.e., 
8 point nodes along (±1,±1,±1) in addition to 6 point 
nodes on the o, b and c axes, which do not agree with the 
cxpcrinicntal data when the magnetic field is tilted from 
the a6- plane. Furthermore, since d\ has all dx, dy 
and dz components in the strong spin-orbit coupling case 
for Oh symmetry, generally di ■ dm 7^ for the second 
component dm, and it is difficult to satisfy constraint 
(i). Only when dm oc C(fe) x d* with arbitrary vector 
function C{k), we can satisfy di ■ dm = 0. However, 
even when the second component dm appears, all point 
nodes for di do not vanish in the /-wave cases, which is 
contrary to that suggested by the experimental data.^' 
Thus, there is no suitable state in this category. 

In the last category of the weak spin-orbit coupling 
triplet state, we can freely select the spin part. How- 
ever, other than the pairing functions already listed in 
the strong spin-orbit coupling case, there are no suit- 
able functions which have point nodes in the Oh symme- 
try among the list of the classifications. However, if we 
consider the ad hoc state, we can obtain some states to 
simultaneously satisfy (i)-(iii) such as 



di(fe) = x{ka ± ikb){kb ± ikc){kc ± ika), 



(24) 



which has 6 point nodes on the a, b, c axes. It can be 
generally concluded that we have to take a broken time- 
reversal symmetry state (i.e., complex function) to have 
point nodes in triplet pairing. 

Having found possible states for the H-phase, such as 
in eq. (24), we now discuss the L-phase, which should 
reduce the symmetry from fourfold to twofold in the ab- 
plane. Below the second transition, the pairing function 
dm satisfying di ■ dm ~ starts to grow as a second- 
order transition, which breaks the fourfold-symmetry in 
the a6-plane. This can be possible if we choose, for ex- 
ample, dm{k) ~ yka or dm{k) — y{kaiikc). The 4 point 
nodes on the a and b axes in the H-phase changes into 
2 point nodes only on the 6-axis, reducing to twofold- 
symmetry in the L-phase. It is noted that the L-phase 
is necessarily non-unitary, where time-reversal symme- 
try is broken in both the orbital and spin parts of the 
pairing function, because d x d* ^ for d = di + ^2, 
as di is a complex function. In this ease, one branch of 
the energy gap |dp ± |d x d*| of the non-unitary state 
has a small gap, but it remains finite because d^ 7^ dy in 
this ease. Therefore, the half-residual density of states 
in the non-unitary states vanishes at low temperature. 
In this L-phase, we expect anisotropic spin susceptibility 
because the spin component is not isotropic. In the non- 
unitary state, we expect the spontaneous moment in the 
spin part of the pairing state, which can be a key point 
to determining the pairing function. 

In summary, for the fourfold (H-phase)-twofold (L- 
phase) transition in PrOs4Sbi2, we have examined the 
mixing effect due to gradient coupling between multi- 
components in unconventional superconductors by cal- 
culating Hc2 near Tc- We then narrow down the pos- 



sible pairing functions in both the H- and L-phases of 
PrOs4Sbi2 in light of the existing experimental informa- 
tion. If we require that the constraints (i)-(iii) mentioned 

above must be fulfilled simultaneously, the candidates for 
the H-phase are the pairing states, such as those in eq. 
(24). The main points of our results are: (A) Any spin- 
singlet states are excluded due to gradient coupling, if 
the second transition is of the second order. The weak 
spin-orbit coupling triplet pairing state is plausible as a 
pairing function. (B) In the H-phase, it is most likely a 
broken time-reversal symmetry state in the orbital part, 
because of point nodes in triplet pairing. In addition, 
the L-phase becomes necessarily a non-unitary state by 
the combination of two order parameters. 

The authors would like to thank Y. Aoki, K. Izawa 
and Y. Matsuda for helpful discussion. 
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